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Abstract: An approximate perturbed direct homotopy reduction method is proposed and 

Ph ' applied to two perturbed modified Korteweg-de Vries (mKdV) equations with fourth order 

! 

dispersion and second order dissipation. The similarity reduction equations are derived 
to arbitrary orders. The method is valid not only for single soliton solution but also for 
the Painleve II waves and periodic waves expressed by Jacobi elliptic functions for both 

> 

ff^ ' fourth order dispersion and second order dissipation. The method is valid also for strong 

lO ■ perturbations. 
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It is very difficult to study nonlinear phenomena lies in the fact that there are various non- 
linear systems which are usually nonintegrable. For some types of idea cases so-called integrable 
models one may use some types of powerful methods (such as the symmetry reduction method [1], 
the Darboux transformation [2] , the nonlinearization [3] or symmetry constraint method [4] etc) 
to find some kinds of exact solutions thanks to there usually exist infinitely many symmetries. 
However, for real nonintegrable physical systems, there are only a little of symmetries or even 
there is no symmetry at all. In many cases, the nonintegrable sector of a physical system may 
company with some small parameters. In these cases, one may use the perturbation theory to 
treat the problems via different approaches. Among these approaches, the approximate sym- 
metry reduction method may be one of the best ways [5,6]. To find symmetry reductions, one 
may use the classical, nonclassical approaches and/or the Clarkson-Kruskal's (CK's) direct 
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method [7, 8] . The CK's direct method is simplest one and can be used to find many group 
invariant solutions without using group theory. Furthermore, in more general cases, the pertur- 
bations may not be weak at all. For strong perturbations, some other types of approaches, such 
as the homotomy analysis method (HAM) [9] and the Linear [10] and nonlinear [11] nonsensitive 
homotopy approaches etc., have to be used. 

In this letter, we try to combined the CK's direct symmetry reduction method and HAM to 
an approximate homotopy direct reduction approach (AHDRA). 

The celebrated modified Korteweg-de Vries (mKdV) equation appears in many branches of 
nonlinear science. As one form of approximation, the singularly perturbed form, 

ut + Gau^Ux + Ur^s = e{u^2 + u^i), (1) 

where a = 1 or a = —1, the subscripts mean the differentiations with respect to x in n times, 
has arisen in a number of physical fields, such as models of shallow water on tilted planes [12]. 
Soliton perturbation property of the mKdV equation was analyzed in [13-15]. 
In this letter, we consider two special forms of the above equation 

Ut + Qau^Ux + Uxxx = eux3±i , (2) 

with fourth order dispersion (the up sign case) and second order dissipation (the lower sign 
case) in terms of APDRA which is a combination of perturbation theory, direct method and 
HAM [7-9]. 

It should be emphasize that in real physical case, the perturbation terms, say, the parameter 
e in ([2]), may not be small. When the perturbations are not weak, the HAM may be successfully 
applied by introducing a homotopy H{u, q) = of the original model A{u) = 0. When the 
homotopy parameter, q = 0, the homotopy model Hq(u) = H(u,0) = should be solved via 
known approach. Usually, Hq{u) is selected as a linear system. In this paper, we select Hq{u) as 
an integrable nonlinear system. Concretely, for the perturbed mKdV system ([2]), we introduce 
the following linear homotopy model (linear for the homotopy parameter), 

(1 - q){ut + 6au^Ux + Uxxx) - q{ut + Gau^u^ + u^xx - (^u^3±i) = 0. (3) 

It is clear that when g = 0, ([3]) is the well known integrable mKdV equation which can be solved 
via many methods. When q = 1, it is just the original model ([2]). Now we can solve ([3]) via 
perturbation approaches by taking g as a perturbation parameter no matter e is small or not. 
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For Eq. ([3]), according to perturbation theory, the solution can be expressed as a series of q 

oo 

u = ^q^Uj, (4) 
j=0 

with Uj being functions of x and t. Substituting Eq. ^ into Eq. ([2]) and vanishing the 
coefficients of ah different powers of ^, we get 

0(e°) : uot + Qau^UQx + Uq^z = 0, (5a) 
0{e^) : uit + (Sa{u^uix + 2uqUiUqx) + ^ix^ ~ eno^^sii = 0, (5b) 
O(e^) : U2t + Qa{uQU2x + u\uqx + '2uqU2Uqx + '2uqUiU2x) + U2x3 - £""1x3*1 = 0, (5c) 

j k 

0{e^) : Ujt + 6a ^ ^ uiUk-iUj^k,x + UjxS - eUj_i^x'-i±i = 0, (5d) 

k=0 1=0 

The similarity solutions for the above equation are of the form 

Uj = Ujix,t,Pjiz{x,t))), (j = 0,l,---), (6) 

where Uj, Pj and z are functions with respect to the indicated variables and Pj{z) satisfy a 
system of ordinary differential equations, which can be obtained by substituting Eq. ^ into 
Eq. (I5|). After the substitution, it is easily seen that the coefficients for Pj^zzz and Pj^zzPj,z 
are Uj^p.z^ and 2>aUj^p.p^z^ respectively. We reserve uppercase Greek letters for undetermined 
functions of z from now on. Because the functions Pj{z) dependent only on the variable z, then 
the ratios of the coefficients are only functions of z, namely, 

^aUj^P^p/x = ^i,P,^'r,(z), (j = 0, 1, • • • ), (7) 

with the solution Uj = Fj{x,t) + Gj{x,t)e^^^^^'^^^ — > Fj{x,t) + Gj{x,t)Pj{z). Hence, it is 
sufficient to seek similarity reductions of Eq. ([5]) in the special form 

Uj = Ujix, t) + (3j{x,t)Pj{z{x,t)), (j = 0,l,---), (8) 

instead of the general form Eq. ([6j). 

Remark: Three freedoms in the determination of aj{x,t), (3j{x,t) and z{x,t) can be notified: 
(i) If aj{x,t) has the form aj{x,t) = a'j{x,t) + Pj^l{z), then one can take 0(z) = 0; 
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(ii) If f3j{x,t) has the form (3j{x,t) = f3'j{x^t)Q.{z), then one can take Q.{z) = constant; 
(hi) If z{x,t) is determined by Q.{z) = zo{x,t), where 0,{z) is any invertible function, then 
one can take Q{z) = z. 

Substituting Eq. ([8]) into Eq. (5a), we can see that the coefficients for Pqzzz, PozPqi PqzPq 
and Pqzz are (3oz^, Ga0QZx, 12aao/3o-^x and S^qxz'^ + 3(3oZxZxx, respectively. We require that 

GaP^Zx = pozl^^oiz), Uaaop^Zx = l3ozl<^o{z), SPoxzl + S/SqZxZxx = l3ozlno{z), (9) 

and thus, applying Remark (i), (ii) and (iii), we have 

ao{x, t) = 0, /3o(x, t) = Zx, z{x, t) = e{t)x + a{t). (10) 

Eq. (5a) is then simplified to 

e^Pozzz + GO^P^Poz + {Otz - eta + eat)Poz + OtPo = 0. (11) 

From the coefficients of Pqzzz-, zPqz and Pqz of the above equation, it is easily seen that 

9t = Ae\ -9ta + eat = Be\ (12) 

with A and B being arbitrary constants. 
When ^ / 0, Eq. has the solution 

1 B 1 

e = -{3A{t - to))~3, = + soit - to)~3, (13) 

where sq and to are arbitrary constants. 

On substitution of Eq. ([8]) into the general form Eq. (5d), the coefficients of Pj^z'h -Pj-i.^sii 
and PqzPq are PjZx, ^ePj-iZx and l2aajPQZx, respectively, which lead to 

Pj-izt = Te^jzl^jiz), 12aajP^Zx = Pjzl^j{z), {j = 1, 2, • • • )■ (14) 

By Po = Zx, Remark (i) and (ii), we obtain 

a, =0, /3j = zl^^, (j = 0,l,---). (15) 

Eq. (HI), (HI), ([To]) . (fT3|) and (fT5]) determine the perturbation series solution to Eq. ^ 

oo 

u = 5;(-l)^+i(3^(t - to))-s^'^^h^P,{z), (16) 

j=0 
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with the similarity variable z = —x{3A{t — to)) +so{t — to) ^ — ^ and the similarity reduction 
equations are 



j k 



-{j + l)APj - {Az + B)Pj^, - 6a 5^ PiPk^iPj-k^z + ei^j-i,.3±i, (j = 0, 1, • • • ), (17) 

fc=0 1=0 

with P_i = 0. When j = 0, Eq. p!7|) degenerates to Painleve II type equation. 
When ^ = 0, Eq. ^ has the solution 

e = tQ, a = B4t + so, (18) 

where sq and to ^ire arbitrary constants. From Eq. (llOp . Eq. (|18p implies an equivalent travelling 
wave form z = x + ct, so that we obtain the perturbation series travelling wave solution to Eq. 

m 

oo 

" ' z = x + ct, (19) 



u = Y.e^Pj{z), 

j=0 



where all Pj{z) satisfy 



j k 

cPj + 2aY,Yl Pj'kPk-lPl + Pj,z^ - ^Pj-i,z2±i + aj =0, (j = 0, 1, • • • ), (20) 

A:=0 /=0 

with P_i = 0. 

Taking j = 0, it is obvious that Eq. (|2U|) becomes the zeroth order equation, the well known 
mKdV equation which has the general solution 

"^0 dp 



±(z - Zo) 



(21) 



V^co - ap"^ - cp2 - 2aop 

with arbitrary constants c, oq, cq and zq. It is also interesting that for the series travelling wave 
solution, it is not difficult to find the solution of Pj can be expressed by Po> result reads 

dp 



P. 



Co - 2aoPo - cPq - aiif 



\J (co — ap''^ — — 2aop')^ 
where cij and C2j are arbitrary constants while 



(co - ap4 - cp2 - 2aop)3 
dpdp' , 



(22) 



j-l A; i-1 

F,(Po) = 2a ^ ^ Pj-kPk-iPi + 2aPo ^ - eP,_i,,3±i + a,-. 

/c=l Z=0 /=1 



(23) 



The general solution (j2ip can be rewritten as some types of Jacobi elliptic functions [16] . For 
some special selections of the constants, it can be written as some types of soliton solutions or 
periodic wave solutions, for instance, if we select a = —1, c = ao = 0, cq = A;^ for the up 
sign (the dissipative case), then we have the hyperbolic tangent shape kink soliton solution 

Pq = kiim\i{kx + 2k^t) = kT. (24) 



Remark: The convergence of infinite series solution Eq. (jl6p is superior to the fourth order 
dispersion case (the up sign case), because the general terms of Eq. (|16p become infinitesimal 
for sufficiently large time 

\2,A{t-t^)\ » 1. 

For the infinite series solution (|16p with the lower sign (the dissipative case), the series will be 
convergent for not very large time, i.e. for 

|3^(t-to)| « 1- 

More specifically, for the dark solitary wave solution (|26p. we can easily find the closed forms 
for the higher order homotopy perturbation solutions. Here is a explicit form of the eighth order 
approximate solution (the end point condition (7 = 1 has been used) 

,^ e S^T-2T y 2T + BS^(BT + 1) , 12T + 3B^S'^ - 23(3"^ - 3)S'^ . 

u = kT -\ \ ; -e^ — e'' 

6 48A; 2034A;3 331776A:^ 

3B^S^{BT - 1) - BS^jTB^ + 15 - 3BT - 2B^) - SOT § 
+ 15925248F e+---, (25) 



where 



S = Vl -T2, S = In i- 



1 + r 

It should be emphasized that homotopy approximation convergence quite well not only for weak 
perturbation (small e) but also for strong perturbations. Fig. 1 shows the schematic plots of 
the first five approximants with respect to the orders 1, 2, 4, 6 and 8 respectively from upper to 
lower of the right side of the figure while the parameters are fixed as 

e = 2.9, k = l. (26) 

From the figure, we find that the lines of the sixth order and the eighth order are almost stuck 
together though the "perturbed parameter" e = 2.9 which is not a small one! 



Figure 1: The plots of the perturbed kink sohtary wave solution for the orders 1, 2, 4, 6 and 8 
respectively from upper to lower of the right side of the figure while the parameters are fixed as 
e = 2.9, k = l. 

Similar to the HAM, the APDRA is applicable to other perturbed nonlinear partial dif- 
ferential equations with and without small parameters and it is thought-provoking to explore a 
general principle for the perturbed nonlinear partial differential equations holding similar results. 
Different from the HAM, we take the zeroth order as an nonlinear integrablc system instead of 
a linear one, which largely modified the convergence rate. Here we take the direct method as a 
tool to find the approximate symmetries and symmetry reductions. The similar results can also 
be obtained via approximate classical and nonclassical symmetry reduction approaches which 
may be used to the KdV-Burgers equation [6], the perturbed nonlinear Schrodinger systems [?] 
and the perturbed Boussinesq system [18]. 
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